Enhanced Conductance Through Side-Coupled Double Quantum Dots 



R. Zitko^ and J. Bonca^'^ 
^J. Stefan Institute, SI- 1000 Ljubljana, and ^Department of Physics, 
FMF, University of Ljubljana, SI-1000 Ljubljana, Slovenia 
(Dated: February 2, 2008) 

Conductance, on-site and inter-site charge fluctuations and spin correlations in the system of 
two side-coupled quantum dots are calculated using the Wilson's numerical renormalization group 
(NRG) technique. We also show spectral density calculated using the density-matrix NRG, which 
for some parameter ranges remedies inconsistencies of the conventional approach. By changing 
the gate voltage and the inter-dot tunneling rate, the system can be tuned to a non-conducting 
spin-singlet state, the usual Kondo regime with odd number of electrons occupying the dots, the 
two-stage Kondo regime with two electrons, or a valence-fluctuating state associated with a Fano 
resonance. Analytical expressions for the width of the Kondo regime and the Kondo temperature 
are given. We also study the effect of unequal gate voltages and the stability of the two-stage Kondo 
effect with respect to such perturbations. 
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I. INTRODUCTION 

The advances in micro-fabrication have enabled studies 
of transport through single as well as coupled quantum 
dots, where at very low temperatures Kondo physics and 
magnetic interactions play an important role. A double- 
dot system represents the simplest possible generalization 
of a single-dot system which has been extensively stud- 
ied in the past. Recent experiments demonstrate that 
an extraordinary control over the physical properties of 
double dots can be achieve d^'^i^i'^ , which enables direct 
experimental investigations of the competition between 
the Kondo effect and the exchange interaction between 
localized moments on the dots. One manifestation of this 
competition is a two stage Kondo effect that has recently 
been predicted in multilevel quantum dot systems with 
explicit exchange interaction coupled to one or two con- 
duction channek^'^. Experimentally, it manifests itself as 
a sharp drop in the conductance vs. gate voltage G{Vg)^ 
or as non-monotonic dependence of the differential con- 
ductance vs. drain-source voltage dl / dVds{Vds)^- 

Fano resonances, which occur due to interference when 
a discrete level is weakly coupled to a continuous band, 
were recently observed in experiments on rings with em- 
bedded quantum dotaS and quantum wires with side- 
coupled dotsiS. The interplay between Fano and Kondo 
resonance was investigated using equation of motioniii^ 
and slave boson techniques''^^. 

In this work we study a double quantum dot (DQD) 
in a side-coupled configuration (Fig. QJ, connected 
to a single conduction-electron channel. Systems of 
this type were studied previously using non-crossing 
approximation'^^, embedding techniquc^^ and slave- 
boson mean field theory^'^'^®. Numerical renormal- 
ization group (NRG) calculations were also performed 
recentlyi^, where only narrow regimes of enhanced con- 
ductance were found at low temperatures. We will show 
that when the intra-dot overlap is large, wide regimes of 
enhanced conductance as a function of gate- voltage exist 



at low temperatures due to the Kondo effect, separated 
by regimes where localized spins on DQD are antiferro- 
magnetically (AFM) coupled. Kondo temperatures Tk 
follow a prediction based on Schrieffer- Wolff transforma- 
tion and poor-man's scaling. In the limit when the dot a 
is only weakly coupled, the system enters the "two stage" 
Kondo regimc"'^^'^^, where we again find a wide regime 
of enhanced conductivity under the condition that the 
high- and the low- Kondo temperatures {Tk and re- 
spectively) are well separated and the temperature of the 
system T is in the interval T]^ < T < Tr-. 




FIG. 1: Side-coupled conflguration of quantum dots 



II. MODEL AND METHOD 

The Hamiltonian that we study reads 

H = 5d{nd - 1) + 5a{na - 1) - ^ {d^a^ + a\dc) 

a 

+ -^y^M-V) + — (n^ - 1) 

ha k<j 

(1) 

where Ud — X^o-'^ct'^c ^^"^ "a ~ X^o- '^o-'^o"- Operators dj. 
and a\ are creation operators for an electron with spin 
CT on site d or a. On-site energies of the dots are defined 
by e = (5 — J7/2. For simplicity, we choose the on-site 
energies and Coulomb interactions to be equal on both 
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dots, 5a = 5d = 5 and Ua = Ud = U. Coupling between 
the dots is described by the inter-dot tunnel coupling t^. 
Dot d couples to both leads with equal hopping t' . As 
it couples only to symmetric combinations of the states 
from the left and the right lead, we have used a unitary 
transformation^^ to describe the system as a variety of 
the single-channel, two-impurity Anderson model. Oper- 
ator c[.^ creates a conduction band electron with momen- 
tum k, spin a and energy = —Dcosk, where D — 2t 
is the half-bandwidth. The momentum-dependent hy- 
bridization function is Vd{k) = —{2/\/N + l)t' sin fc, 
where N in the normalization factor is the number of 
conduction band states. 

We use Meir-Wingreen's formula for conductance in 
the case of proportionate coupling^*' which is known 
to apply under very general conditions (for exam- 
ple, the system need not be in a Fermi-liquid ground 
state) with spectral functions obtained using the NRG 
technique2ii22i2^i24. At zero temperature, the conduc- 
tance is 

G = GQT:Tpd{0), (2) 

where Go = 2e^/h, Pd{'-^) is the local density of states of 
electrons on site d and T/D — (t' /t)^. 

The NRG technique consists of logarithmic discretiza- 
tion of the conduction band, mapping onto a one- 
dimensional chain with exponentially decreasing hopping 
constants, and iterative diagonalization of the resulting 
Hamiltonian^"'^. Only low-energy part of the spectrum 
is kept after each iteration step; in our calculations we 
kept 1200 states, not counting spin degeneracies, using 
discretization parameter A — 1.5. 

III. STRONG INTER-DOT COUPLING 

A. Conductance and correlation functions 

In Fig. [2Ji we present conductance through a double 
quantum dot at different values of intra-dot couplings vs. 
d. Due to formation of Kondo correlations, conductance 
is enhanced, reaching the unitary limit in a wide range 
of 5. 

To better understand multidot problems in the case 
of strong inter-dot coupling, it is helpful to exactly di- 
agonalize the part of the Hamiltonian that describes the 
dots and rewrite the entire Hamiltonian (Eq.^ in a form 
similar to the ionic modeiSS.: 

77„ = ^e,4^^Cfc,. + ^i?(a)|a)(a| (3) 

Hi= ik.\p^a \a){l3\ Ck,a+tl^\p^J^, \P){a\. (4) 

where multi-indeces a and (3 stand for quantum num- 
bers {Q, S, Sz,r). Here Q is the charge number Q = 




5/D 



FIG. 2: Conductance and correlation functions of DQD vs. 
S. Besides different values of td, indicated in the figure, other 
parameters of the model are T/D = 0.03 and U/D — 1. Tem- 
perature is T/D = 10~^, which for all parameters used cor- 
responds to a zero temperature limit. In particular, in the 
Kondo plateaux T <C Tk for all S. 



(ua — 1) + {ud — 1), S, Sz are the spin and its component, 
while r numbers different states with the same Q, S, Sz 
quantum numbers. In the absence of the magnetic field, 
^2 is irrelevant and will be omitted from now on. Each 
(Q, S, r) multiplet is then 25-1-1-fold degenerate. Finally, 
the effective hopping coefficients 

tka\p^c.^Vd{k){a\dl\(3) (5) 

correspond to electrons hopping from the conduction 
band to the dots. 

The eigenvalue diagram in Fig. represents the gate- 
voltage dependence of the multiplet energies E{Q, S, r). 
From this diagram we can read off the ground state and 
the excited states for each parameter S. 

Regimes of enhanced conductance appear in the in- 
tervals approximately given by < \S\ < 62 ^ where 
di = td{2^\ + {V/MdY - 1) and<52 = {V /2 + td). These 
estimates are obtained from the lowest energies of states 
with zero, one and two electrons on the isolated double 
quantum dot: 

£;(-2,0,0) [/- 25, 

£;(-l,i,0) = C//2-5-t<i, (6) 
£;(0, 0, 0) = -2td^J\ + (IJ/Mdf + C//2. 

The widths of conductance peaks (measured at G/Gq = 
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measured by An^ 



as shown in Fig.[2li. Fi- 



FIG. 3: Eigenvalue diagram for isolated DQD system. The 
diagram is symmetric, since for S —5, E{Q,S,r) 
E{—Q,S,r). Points A and B correspond to valence- 
fluctuation regions when the charge on the dot changes, while 
point C corresponds to the center of the Kondo regime, 
when the Kondo temperature is the lowest. Parameters are 
U/D = 1 and td/D = 0.3. 



1/2) are therefore approximately given by 



A = U/2 + 2td{l - ^l + iU/Ud)^). 



(7) 



Note that in the Umit of large td, A ~ U/2, and in the 
limit of large [/, A ~ 2td- As it will become apparent 
later in this paper, this naive estimate fails when ^ 0. 
Comparison of conductance-peak widths with the ana- 
lytical estimate A is shown in Fig. 
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FIG. 4: a) The width of conductance peak A vs. td as ob- 
tained from NRG calculations (full circles), compared with 
the analytical result as given in the text (dashed line), b) 
Kondo temperatures Tk vs. 5 as obtained from the widths 
of Kondo peaks (full circles). Analytical estimate, Eq. |H1 is 
shown using dashed lines. The rest of parameters are identical 
to those in Fig. |5| 



nally, we show in Fig. |2t spin-spin correlation function 
(Sa ■ Sd). Its value is negative between two separated 
Kondo regimes where conductance approaches zero, i.e. 
for —Si<S<Si, otherwise it is zero. This regime further 
coincides with (n) ^ 2. Each dot thus becomes nearly 
singly occupied and spins on the two dots form a local sin- 
glet due to effective exchange coupling Jeff w 4i^/J7. Sur- 
prisingly, the absolute value of (Sa • S^) increases and it 
nearly reaches its limiting value, i.e. {Sa ■ S^) = —3/4, as 
the intra-dot hopping td, and with it Jeff, decrease. This 
seeming contradiction is resolved by recognizing that Jeff 
represents the exchange interaction based on the effec- 
tive Heisenberg coupling between localized spins on DQD 
only in the limit when td/U — *■ 0. Smaller absolute val- 
ues of (Sa • Sd) are thus due to a small amount of double 
occupancy at larger values of td/U. 

In Fig. we present Kondo temperatures Tk vs. S 
extracted from the widths of Kondo peaks. Numerical re- 
sults in the regime where (n) ^ 1 and 3 fit the analytical 
expression obtained using the Schrieffer- Wolff transfor- 
mation that leads to an effective single S" = 1/2 Kondo 
model. We obtain 



Tk = 0.182f7v/A^cxp[-l/poJ] 



(8) 



with 



2r 

poJ = — 



a 



where 



TT \\E{-l,^,0)-E{-2,0,0)\ 

, i ] 

\E{-l,lO)-E{0,0,0)\J' 

a=|(-l,i,i,0|d||-2,0,0,0)|2 = l/2, 
/3= |(0,0,0,0|4|-l,i,i,0)|2 

Ud + U+ ^Wtj + [/2) 



(9) 



(10) 
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(i6t2 + u[u + y^mtf+TP 



The factor 0.182t/ in Eq. |H1 is the effective bandwidth. 
The same effective bandwidth was used to obtain Tk of 
the Anderson model in the regime U < £J?M6:. 



Next, we focus on various correlation functions. In 
Fig. we show S, calculated from expectation value 
(S^^t) = S{S + 1), where Stot = Sa + S^ is the total 
spin operator. S reaches value 1/2 in the regime where 
G/Gq — 1. Enhanced conductance is thus followed by 
the local moment formation, indicative of the Kondo ef- 
fect. This is further supported by the average double-dot 
occupancy (n), where n = Ua + n-d, which in the regime 
of enhanced conductivity approaches odd-integer values, 
i.e. (n) — 1 and 3 (see Fig. |2t). Transitions between 
regimes of nearly integer occupancies are rather sharp; 
they are visible as regions of enhanced charge fluctuations 



B. Spectral densities 

We observed that the spectral density calculation us- 
ing the conventional NRG approachSi fails for our model. 
The spectral densities manifest spurious discontinuities 
and the normalization sum rule is violated for some 
choices of model parameters (Fig. E)). This happens be- 
cause at the early stages of the NRG iteration the low- 
est energy state does not yet correspond to the true 
ground state. We found that we obtain correct results us- 
ing the density-matrix NRG technique^^, which remedies 
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FIG. 5: Comparison of spectral densities computed by the 
conventional NRG and by the DM-NRG for parameters where 
the conventional NRG approach has difficulties: td/D = 0.3, 
U/D = 1, r/D = 0.03, 5/D = 0.48. Note the discontinu- 
ities in the spectral density calculated using the conventional 
NRG. At low frequencies both method produce identical re- 
sults, so we do not show the Kondo peak which exceeds the 
vertical scale of the plot. In this calculations, the effective 
temperature is zero. 



the shortcomings of the conventional approach. We suc- 
ceeded in implementing this technique with eigenstates 
defined within subspaces with well defined occupation Q 
and total spin S, as opposed to well defined occupation 
and total spin component Sz (see the Appendix). The 
improvement in numerical efficiency is sufficient to en- 
able consideration of more complex systems, such as the 
double quantum dot. 
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FIG. 6: Zero-temperature spectral density pd{i^) sweeps for 
td/D — 0.3. a) Spectral density calculated using the conven- 
tional NRG approach, b) Spectral density calculated using 
the density-matrix NRG approach. Note that the vertical line, 
representing the Kondo peak, has been artificially broadened. 
Its true width is given by Tk- 

In Fig. El we present sweeps of Pd{^) calculated us- 
ing both approaches. In vast regions of the plot the re- 
sults are in perfect agreement. The differences appear for 
those values of S where the ground state changes. Three 
characteristic spectral density curves calculated using the 
DM-NRG are shown in Fig. 

Features in the spectral density sweeps can be easily 
interpreted using eigenvalue diagram in Fig. |21 At low 
temperatures and for constant S, spectral density A(uj) 
will be high whenever the energy difference AE = Ei—Eq 




FIG. 7: Three zero-temperature spectral densities pd{ui) for 
td/ D = 0.3: at the particle- hole symmetric point = 0, in 
the Kondo regime 5/D = 0.62 and in the empty-orbital regime 
5/D — 1. Inset: Scaling of spectral densities pd{u>/TK) in the 
Kondo region 0.52 < S < 0.76. Parameters are as in Fig. |5| 



between the ground state (0) and an excited state (1) is 
equal to +uj (particle excitations, Qi = Qo + 1, Si ~ 5o± 
i) or to — cj (hole excitations, Q\ = Q{y — \, S\ = 5*0 ±i). 
At (5 = two broad peaks are seen located symmetrically 
at w - ±(5i (see Fig. Eland Fig. 13 at 8/D = 0). At this 
point the model is particle-hole symmetric and therefore 
E{Q, S, r) = E{-Q, S, r) for all' Q, S, r. Consequently, 
the spectrum is also symmetric, Pdi^) ~ Prf(~^)- With 
increasing 6, the particle excitation energy E{1, ^,0) — 
i?(0, 0,0) increases and the corresponding peak quickly 
washes out. The hole excitation energy i?(— 1,^,0) — 
E{0, 0, 0) decreases and the peak gains weight. 

At 5 = (5i, £;(0, 0, 0) = E{-1, i, 0) (point A in Fig. O 
and the system enters the Kondo regime: a sharp many- 
body resonance appears which remains pinned at the 
Fermi level throughout the Kondo region (see Fig. El and 
Fig. El at 6/D = 0.62). Kondo effect occurs whenever the 
ground state is a doublet, S — 1/2, and there are excited 
states with S' = 0, 1, Q' — Q ± 1. We recognize such 
regions by characteristic triangular level crossings in the 
eigenvalue diagram, one of which is marked as triangle 
ABC in Fig.O The high-energy peaks at w = E{0, 0, 0)- 
£;(-l,i,0) > and w = E{-1,^,0) - £^(-2,0,0) < 
in the spectral density are also characteristic: they cor- 
respond to particle and hole excitations that are at the 
heart of the Kondo effect. 

In the case of the DQD we also see additional struc- 
ture for 6i < S < 62- a broad peak at w = E{Q, 1, 0) — 
i?(— l,i,0) which corresponds to virtual triplet excita- 
tions from the ground state. These excitations could also 
be taken into account in the calculation of the effective 
exchange interaction, Eq. ||^, however due to their high 
energy, they only lead to an exponentially small renor- 
malization of the Kondo temperature, which may be ne- 
glected. 

In the inset of Fig. [7\ we show scaling of Kondo peaks 
vs. lo/Tk- In the case of perfect scaling, all curves should 
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exactly overlap. However, Kondo temperatures of differ- 
ent peaks differ by almost four orders of magnitudes, as 
seen in Fig. 0)3. Moreover, Kondo peaks become asym- 
metric near the edges of the Kondo region, i.e. for 6 ^ Si 
and S 62- Note also that for each point in Fig.^Ja there 
is a respective spectral density presented in the inset of 
Fig.CI 



C. Effect of the on-site-energy splitting 

In real double quantum dot systems, it is difhcult to 
produce dots and electrodes with identical properties. In 
particular, it is not easy to achieve equal on-site energies 
6d = — S for a wide range of parameters i5. Therefore 
it is necessary to study the robustness of the physical 
regimes with respect to the splitting of the on-site ener- 
gies. 

We now generalize slightly our model to allow for un- 
equal on-site energies by introducing a new parameter, 
K, so that on-site energies are Sd = S + k and 6a ~ 5 ~ k. 
The DQD part of the Hamiltonian is now 

Ho = {S + K){nd -1) + {S- K){na - 1) 



(11) 



l^{nd-lf + -{na~lf 



Only K > has to be considered, because Hq{5, —hi) 
Ho{—S,k) when particle-hole transformation dj^ — s- 
dp, ajj — + — flp is performed. 

By diagonalization Hq exactly, we find that the effect of 
K is to shift states while the general eigenvalue structure 
is maintained and is similar to that shown in Fig. |31 We 
therefore still expect to see two Kondo peaks as S is swept 
for a constant n. 

In the limit of large n and for S such that Sd ^ and 
6a ±00, the side-coupled dot becomes irrelevant and 
we recover the familiar single-impurity Anderson model 
with conduction plateau in the region —U/2 < Sd < U/2 
with Tk given by Eq. © with 



2r 

PoJ = — 



1 



1 



TT \\Sd-U/2\ \Sd + U/2 



(12) 



This scenario is corroborated by NRG calculations 
(Fig. ISJ. The approach to the limit can be determined 
using the expressions Eqs. JHl and ((HJ, with k depen- 
dent a^n), P{k) and Ei{k), E2{k). The result is shown in 
Fig. I^li. The k dependence is non-monotonic, which can 
be traced to non-monotonic behavior of coefficient /?(t). 
At k/D = 1 the asymptotic value Tk/D = 1.04 x 10""^ 
is nearly reached. 

In the limit of large k, but for S such that iJ^ ~ and 
Sd ±00, the embedded dot will be either fully occu- 
pied or empty with very little charge fluctuations, while 
the side-coupled dot will maintain 5 = 1/2 spin moment. 
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FIG. 8; Conductance as a function of gate- voltage Sd for 
different gate-voltage asymmetries k. Label SIAM corre- 
sponds to results for conductance of the corresponding single- 
impurity Anderson model. The temperature is T/D — 10~®. 
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FIG. 9: Kondo temperature as a function of /t for a) the 
Kondo plateau at 5^ ~ and b) the Kondo plateau at Sa ~ 
0. In both cases we calculate the Kondo temperature in the 
middle of the plateau, where Tk is lowest. Parameters are 
U/D = 1, r/D = 0.03 and ta/D = 0.3. 



This moment will be screened, but due to the suppres- 
sion of the virtual electron hopping through the dot d 
the corresponding Kondo temperature rapidly drops with 
increasing k (Fig. Eb)- Keeping the temperature con- 
stant and increasing k, the Kondo plateau in the region 
—U/2 < Sa < u/2 will therefore quickly evolve into two 
peaks separated by a Coulomb blockade valley and for 
very large k it will not conduct at all. This prediction 
can also be verified using NRG fFig. I1(J|) . 
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FIG. 10; Conductance as a function of gate-voltage 5a for dif- 
ferent gate-voltage asymmetries k. The temperature is fixed 
to T/D = 10"'* for all k.. 




FIG. 11: Conductance and correlation functions at td/D — 
0.001 (a,...,e), td/D = 0.01 (f,and g) and a blow-up of f) 
in the inset of g). Dashed lines in a) represent G/Go and 
b) S of a single quantum dot with otherwise identical pa- 
rameters. Dashed line in c) represents (ud) of DQD, and 
finally dashed line in the inset of g) represents the semi- 
analytical model described in the text. In h) a schematic 
plot of different temperatures and interactions is presented 
as explained in the text. NRG values of the gap in pd{u>) 
at LJ = and T <g are presented with open circles and 
squares for td/D = 0.001 and 0.01 respectively. Values of J^s 
and analytical results of are presented with dashed and 
dot-dashed lines for td/D = 0.001 and 0.01 respectively. For 
analytical estimates of different values of a = 2.2 and 1.3, 
respectively, were used. Other parameters of the model are 
r/D = 0.03, C//D = l. 



IV. WEAK INTER-DOT COUPLING 
A. Conductance and correlation functions 

We now turn to the limit when 0. Unless oth- 

erwise specified, we choose the effective temperature T 
to be finite, i.e. T ~ 10~^D, since calculations at much 
lower temperatures would be experimentally irrelevant. 
In this case one naively expects to obtain essentially iden- 
tical conductance as in the single-dot case. As S decreases 
below 5 ~ U/2, G/Gq indeed follows result obtained for 
the single-dot case as shown in Fig. Illb . In the case 
of DQD, however, a sharp Fano resonance appears at 
6 — U/2. This resonance coincides with the sudden jump 
in S, (n), as well as with the spike in An^, as shown in 
Figs lllb .c. and d, respectively. Fano resonance is a con- 
sequence of a sudden charging of the nearly decoupled 
dot a, as its level e crosses the chemical potential of the 
leads, i.e at e = 0. Meanwhile, the electron density on 
the dot d remains a smooth function of 6, as seen from 
(nd) in Figlllb. With increasing td, the width of the reso- 
nance increases, as shown for td/D = 0.01 in Fig. lilt'. For 
td ^ 0.1, the resonance merges with the Kondo plateau 
and disappears (see Fig.^Jt). 

The resonance at 5 ~ U/2 can be explained using a 
simple semi-analytical model. For td — 0, we write the 
exact Green's function of the impurity d aX lu — a,s Gd — 
l/(ir— r cot </))2i, where is the scattering phase shift for 
single impurity model that we calculate using NRG. The 
side-coupled dot is taken into account using perturbation 
theory. We obtain full Green's function from the Dyson's 
equation 

gd^g"d/ii-tig',g"d)- (13) 

For the Fano resonance a,t e — S — U/2 ^ 0, we keep only 
the low-energy pole in the Green's function, 

g» « zj^,, (14) 

where z = ^ for e < and z ~ 1 for e > 0. The conduc- 
tance is 

2T\'sm^<P 

° z2t4 + 2r2e2 - zHj cos(2,^) - 2ztleT sin(20) ' 

(15) 

Results of the NRG calculation are compared to the 
prediction from Eq. ^] in the inset of Fig. Illb . We see 
that general features are adequately described, but there 
are subtle differences due to non-perturbative electron 
correlation effects. Numerically calculated Fano reso- 
nance is wider than the semi-analytical prediction and 
G/Gq does not drop to zero. In particular, from Eg. 1151 
it follows that G = at (5 = C//2 (e = 0) and G = Go at 
,5 = U/2 + titan (j}/T. This is not corroborated by NRG 
results, which show maximal conductance at S ~ U/2. 

We now return to the description of the results pre- 
sented in Fig. Illb in the regime where 6 < U/2. As 5 
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further decreases, the system enters a regime of the two- 
stage Kondo effect^''. This region is defined by Jes < Tx 
(see also Fig.lllh^. where Tk is the Kondo temperature, 
approximately given by the single quantum dot Kondo 
temperature, Eq. |Slwith 

2r / 1 1 \ 

Just below 5 < U/2, T falls in the interval, given by 
T^i<$:T<$:TK, where 

T°, ^ Tk exp{-aTK/Jcs) (17) 

denotes the lower Kondo temperature, corresponding to 
the gap in the spectral density Pd{^) at w = and a is of 
the order of one^^. Note that NRG values of the gap in 
Pd{i^) (open circles and squares), calculated at T <C T]^, 
follow analytical results for T^{6) when Jog < Tr-, see 
Fig. lllh . while in the opposite regime, i.e. for Jeff > Tk, 
they approach Jeff. 

As shown in Fig. ^ for 0.3D <S< U/2, G/Gq cal- 
culated at T = IQ^^D follows results obtained in the sin- 
gle quantum dot case and approaches value 1. The spin 
quantum number S in Fig. lllb reaches the value S ^ 0.8, 
consistent with the result obtained for a system of two 
decoupled spin- 1/2 particles, where (S^) = 3/2. This re- 
sult is also in agreement with (n) ~ 2 and the small value 
of the spin-spin correlation function (Sq • S^), presented 
in Fig.im andlTTfe respectively. 

With further decreasing of 5, G/Gq suddenly drops to 
zero ai 6 ^ 0.3D. This sudden drop is approximately 
given by T < T^{d), see Figs. Hit and h. At this point 
the Kondo hole opens in Pdi^) at a; = 0, which in turn 
leads to a drop in the conductivity. The position of this 
sudden drop in terms of S is rather insensitive to the 
chosen T, as apparent from Fig. Illh . 

Below S < 0.251?, which corresponds to the condi- 
tion Jeff ^ Tk{S), also presented in Fig. [TTh . the sys- 
tem crosses over from the two stage Kondo regime to a 
regime where spins on DQD form a singlet. In this case 
S decreases and (Sq • S^) shows strong antiferromagnetic 
correlations. Figs. Illb and e. The lowest energy scale 
in the system is Jeff, which is supported by the obser- 
vation that the size of the gap in Pd{^) (open circles in 
Figs. [TTh ^ is approximately given by Jeff. The main dif- 
ference between td/D = 0.001 and td/D = 0.01 comes 
from different values of Jeff = 4t^/C7. Since in the lat- 
ter case Jeff is larger, the system enters the AFM singlet 
regime at much larger values of (5, as can be seen from 
comparison of Figs.ITT^ and f. Consequently, the regime 
of enhanced conductance shrinks. 



B. Effect of the on-site-energy splitting 

We again explore the effect of unequal on-site energies. 
Unlike in section Fill CI we now introduce the asymmetry 
so that only the level a is shifted by parameter k: 5d = 5 



and 5a = 5+k. We are mainly interested in the parameter 
range where for k — Q the system exhibits the two-stage 
Kondo effect. 

This study will be conducted by computing the impu- 
rity susceptibility 

ximp(r) = {{SD - {sl),) , (18) 

where the first expectation value refers to the system with 
the double quantum dot, while the second refers to the 
system without the dots. In traces of Tximp, the two- 
stage Kondo effect manifests as two successive decreases 
of the susceptibility, at T ^ Tk from around 0.5 to 0.25 
and at T ~ from around 0.25 to 0^^ We now inves- 
tigate how this behavior changes when k 7^ by compar- 
ing temperature dependent susceptibilities calculated for 
a range of parameters At. 




T/D 



FIG. 12; Dynamic magnetic susceptibility ksTxiT) /(gps)^ 
for different k. Parameters are U/ D — 1,T / D — 0.03, td/D = 
8 X 10~^ and 5/D = 0. These calculations were performed 
using a different NRG discretization scheme— with A = 4 and 
keeping 600 states. 

We find that the two stage Kondo effect is robust 
against k ^ perturbation. In fact, it survives until 
K = U/2, see Fig. ^2 For k < U/2 the main effect of 
K is to reduce the lower Kondo temperature T]^. For 
K > U/2, dot a becomes irrelevant at low temperatures 
and dot d is in the usual Kondo regime. 

The transition from k < U/2 to k > U/2 is continuous 
despite a rather abrupt change in Ximp caused by a rapid 
change in occupation of dot a that occurs with increasing 
K on an energy scale proportional to t^. For k = U/2, 
£a = 5a — U/2 = and the side-coupled dot is in the va- 
lence fluctuation regime. For two uncoupled dots, one in 
local-moment regime and the other in valence-fluctuation 
regime, the magnetic moment ksTx/ {gpsf' is expected 
to be 1/4-1-1/6 = 5/12. At high temperatures, x is indeed 
in the vicinity of this value. 
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V. CONCLUSIONS 

In this paper we have explored different regimes of the 
side-coupled DQD. When quantum dots are strongly cou- 
pled, wide regions of enhanced, nearly unitary conduc- 
tance exist due to the underlying Kondo physics. An- 
alytical estimates for their positions, widths, as well as 
for the corresponding Kondo temperatures, are given and 
numerically verified. When two electrons occupy DQD, 
conductance is zero due to formation of the spin-singlet 
state which is effectively decoupled from the leads. In this 
regime most physical properties of the coupled DQD can 
be predicted using the eigenvalue diagram of the isolated 
DQD system. This happens because energy separations 
of different DQD states are much larger than the Kondo 
temperature. Introducing the on-site-energy splitting be- 
tween DQD induces new, even though not unexpected 
effects within the Kondo regime. With increasing level- 
splitting, the Kondo plateau is expanded when the level 
of the embedded dot is kept near the local particle-hole 
limit. In the other case, when, in contrast, the level of 
the side-coupled dot is kept near the particle-hole limit, 
Kondo plateau splits at finite T due to abrupt decrease 
of Tk with increasing level-splitting. 

When quantum dots are weakly coupled Fano resonance 
appears in the valence fluctuation regime. Its width is en- 
hanced as a consequence of interactions which should fa- 
cilitate experimental observation. Unitary conductance 
exists when two electrons occupy DQD due to a two- 
stage Kondo effect as long as the temperature of the 
system is well below Tk and above T^. The experi- 
mental signature of the two-stage Kondo effect in weakly 
coupled regime should materialize through the inter-dot- 
coupling sensitive width of the enhanced conductance vs. 
gate voltage. Experimental observation of the two-stage 
Kondo effect would require no external magnetic field as 
used in Ref/. The inter-dot coupling strength can be 
experimentally varied using gate electrodesi*2i2i4. 
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APPENDIX: DENSITY-MATRIX NRG IN (Q,S) 
SUBSPACE BASIS 

Density-matrix numerical renormalization group (DM- 
NRG)24 was originally implemented for a NRG iteration 
that uses subspaces with well defined charge Q and spin 
projection Sz and it was used to determine the effect of 
magnetic field on spin-projected spectral densities. In 
the absence of magnetic field, when rotational invariance 
in the spin space is recovered, it is advantageous to use 
subspaces with well defined charge Q and spin S, which 
leads to greatly enhanced numerical efficiency. 

Here we show how DM-NRG can be implemented in 



(Q,S) subspace basis in the case of a single conduction 
channel. 

Density-matrix at the last iteration step (estimated us- 
ing the truncated basis) is: 

P^^ exp{-(3EQSu.)\QSSzU;){QSSM, (A.l) 

where cu enumerates different states in each (Q, S) sub- 
space and the grand-canonical statistical sum Z is 

Z = Tr [exp(-/3i7)] = ^ exp(-/3^Q5c.). (A.2) 

Unitary transformation of states from to + 1-th 
stage is 

\QSSzUj)n+i = ^^UQs{iLi\ra)\QSSzra)N+i, (A.3) 



where Uqs are the corresponding unitary transformation 
matrices obtained during the diagonalization step of the 
NRG and the \QSSzCtr)N+i states are defined as^i 

\QSSzrl)N+i = \Q + l,SSzr)N 

\QSSzr2)N+i = uflj^^^^^\Q,S - -,Sz - ^,r)N 

\QSSzr3)N+i = S+^,Sz- ^,r)jv 

+ yflN+i)i\Q'^ + ^,5"^ -I- 5,r)Ar 

\QSSzrA)N+i = /(V+iji^V+ijilQ ^ l,SSzr)N\ Ti), 

(A.4) 

where fj^jy^i-^^ is the creation operator for electrons on 
the {N + l)-th site of the "hopping Hamiltonian" and 
M, V, w, y are the Glebsch-Gordan coefficients 



S + Sz 
25 

S — Sz 
2S 



1/2 



1/2 



S-Sz + l 
25 + 2 



1/2 



y 



S + Sz + l 



1/2 



25 + 2 

The density matrix in the basis of \QSSzar) Ar-|_i states 

is 

P= ^^v{~I3Eqsu) Y UQs{i^\ra)UQs{i^\r' a') 

QSSz<^ ra,r'o(' 

\QSSzra){QSSzra\. 

(A.5) 
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We now perform a partial trace over the states on the additional {N + l)-th site to obtain projector operators 
defined on the chain of length N. Diagonal projectors (those with a = a') are: 

TiN+i {\QSS,rl){QSSyi\) = \Q + 1, SS,r)N{Q + 1, SS^In 
TiN+i{\QSS,r4){QSS,r'4\) = \Q - 1, SS,r)N{Q - I, SS,r\N 



TiN+i i\QSS,r2){QSSy2\) = u^\Q, 5 - i, 5, - i, r)N{Q, 5 - i, 5, - i, rV 

+ v^\Q,S-^,S, + ^,r)N{Q,S-^,S, + ^,r'\N 

Ttn+i {\QSS,r3){QSS,r3\) =w^\Q,S + ^,S, - i,r);v((3,5+ - i,rV 

+ y^\Q,S+lS, + lr)N{Q,S+l,S, + lr'\N 

The out-of-diagonal a = 2, a' = 3 terms give zero when summed over 

s 



(A.6) 



^ TiN+i\QSS,r2){QSS,r'3\ = 



-s 

s 



^ ^ (_ l {S + S,){S-S, + l) I {S - + 1){S + S,) \ 
s^ir^+iV V 25(25 + 2) +V 25(25 + 2) ) 

10, 5 - i, 5, - i, r)w(g, 5 + i, 5, - i, r'U = 

where we took into account that vss = and us,-s = in order to expand the range of summation index Sz- 
Similarly we show that a = 3, a' = 2 terms drop. 

The a = a' = 2 terms can be simplified after summing over 5^: 



^ Trjv+i(|055,r2)(055,r'2|) 

Sz=-S 

s 



Sz=-S+1 

Sz=-(s-h 



(A.8) 



N 



The spin multiplicity of {QS)n+i space is 25+1, while 
the spin multiplicity of {Q,S — ^)Ar is 25. The factor 
(25+ l)/(25) is therefore merely a normalization factor. 
In the last line we emphasized that in the A''-site space 
the Sz runs over all permissible values for spin 5 — ^ . 

By analogy we show that the a = a' = 3 terms also 
simplify: 

s 

^ Tr^+i(|Q55,r3)(g55,r3|) 



Non-zero partial traces of projector operators are 
therefore: 



^Trw+i|055,ra)(055 

Sz 



(A.10) 



JV 



25+1 
25 + 2 



J2 \Q,S+^,Sz,r)M{Q,S+^,Sz,r' 



Sz=-{S+^) 



N 



(A.9) 



Again, the spin multiplicity of the ((5,5+ |)/^ space 
is 25 + 2 and the factor (25 + l)/(25 + 2) takes care of 
the correct normalization. 



with ci 

1, Q2 = 



= C4 



■ I, C2 -- 

Q, Qi 



2S+1 



2S 

Q- 



C3 = 
1, Si 



2S+1 
25+2' 
= 54 



Qi 



Q + 

S, S2 = 

S — ^, S3 = S + ^ and corresponding 52^ ranges over 

all possible values for a given 5q. The reduced density 
matrix remains diagonal in its {QS) subspace index. 
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In general we therefore have 

QSS^ ujLu' 

QS Louj' ra.r' 
UQs{i^'\r'a)Ca{S) E \QaSaS:,ar) iq 



N- 

(A.ll) 



This is to be comapred with 



Pr'^duccd = E E C^r'-''\QSS.r){QSSy\. (A.12) 

QSS^ rr' 

We finally obtain the recursion relation for calculation 
of coefficients C5?,f in the reduced density matrix: 



UJUJ' 

+ E Cl+^-^'^+ VQ+i,5(c.|r4)C/Q+i,s(c.'|r'4) 



This is the main result of the derivation. Using known 
Uq^s matrices, recursion in Eq. (| A. 13(1 is applied after 
the first NRG run to calculate reduced density matrices 
for all chain lengths. In another NRG run, the spectral 
density functions are then calculated with respect to the 
reduced density matrices: 



reduced 
im 



+ (z|4|m)0l4|m)p5f-<=^)<5 (c. - {E, E^)) 

(A. 14) 



2S + 2 Q,S+^,N+l 
25' + 1 ^ t^o;' 

UJUj' 

C/^ i(c^|r2)C/^ i(a;'|r'2) 



Q,S+2 

2S 



2S 



2^ / y i^uj' 



S+7 
,S-\,N + l 



(A.131 
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